ON THE CHOQUET-BRUHAT YORK FRIEDRICH FORMULATION OF THE 

EINSTEIN-EULER EQUATIONS 
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Cn . Abstract. Short-time existence for the Einstein-Euler and the vacuum Einstein equations is proven 

^h ' using a Friedrich inspired formulation due to Choquet-Bruhat and York, where the system is cast into 

C*H, a symmetric hyperbolic form and the Riemann tensor is treated as one of the fundamental unknowns 

of the problem. The reduced system of Choquet-Bruhat and York, along with the preservation of the 

gauge, is shown to imply the full Einstein equations. 
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1. Introduction 



Oh 
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In the vast amount of literature that exists on the Cauchy problem of General relativity (GR) 1 , 

the formulation in terms of a first order symmetric hyperbolic system (FOSH) has recently attracted 

significant attention (see e.g. [YB, Fri2, Fril, FN, FR] and references therein.). Here, we focus on the 

Choquet-Bruhat and York [YB] formulation of the Einstein-Euler system in terms of the Lagrangian 2 

description of the fluid flow, which itself was adapted from an earlier formulation by Friedrich [Fri2]. 

£> \ In [YB] the authors wrote a system of equations in terms of the Riemann tensor (as opposed to the 

ON ' Weyl tensor used in [Fri2]), and chose a gauge that reduced this system to a FOSH. It turned out 

that this system had physical characteristics (in contrast with Friedrich's one), i.e. the assumption 

I/-) \ that the speed of sound in the fluid is less than that of light was crucial to prove the hyperbolicity 

t-h | of the equations. This is important, for example, because it gives a natural breakdown criterion for 

the problem of long-time existence. However, the task of proving that the gauge is preserved and the 

original Einstein-Euler system is satisfied was not carried out in [YB]. In this article, we complete 

the proof of short-time existence for the Einstein-Euler system a la [YB]. In what follows, we shall 

restrict ourselves to barotropic fluids. The more general situation will be addressed elsewhere. 

£j ■ 2. Summary of results 

In the study of the Cauchy problem in GR, one is usually given a Riemannian smooth 3-fold 
(£, /io)j a symmetric 2-tensor K, and other initial data corresponding to the matter fields. This 
initial data is required to satisfy certain constraint equations (essentially the Einstein equations at 
t = combined with the Gauss-Codazzi-Mainardi equations). This prescription of data is usually 
facilitated by means of the conformal method of solving the constraints. The aim then is to find an 
Einsteinian development, i.e. a Lorentzian 4-fold (M,g) = (Ex [0,T],g) containing matter fields 



1 A complete or extensive account of all references is beyond the scope of this paper, whose length we tried to keep 
short. We refer the interested reader to the monographs [B, KN] and the survey papers [CGP, FR]. A long, although 
far from complete, review of the literature of the Cauchy problem for the Einstein-Euler system specifically, is given in 
[D]. 

intuitively, the picture is like this. We can think of two ways to study the flow of a river: one could float downstream 
on a boat, or one could sit on the bank and observe the flow. The former (Lagrangian description) corresponds to 
tracking the position of every particle and the latter (Eulerian description) to observing the velocity vector field. Both 
descriptions are useful in the study of relativistic and non-relativistic fluids. 

1 
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satisfying Einstein's equations with (E, ho) embedded isometrically with K as its second fundamental 
form and obeying the initial conditions on the matter fields. 

Just as in [D], we assume that we are given the aforementioned type of initial data. This data is 
then converted into the type we need for solving the FOSH. In what follows, for the fluid case, p > 
indicates pressure, p{p) > is the density as a function of p, v is the initial 3-velocity of the fluid 
on E, and u is the 4- velocity field on M. We prove our results in the uniformly local Sobolev spaces 
H^ t of Kato [B]. In what follows, E, ho, and p denote the quantities just described, and repeated 
indices are summed over. We also assume the reader is familiar with the terminology and the Cauchy 
problem in GR and the initial conditions for the Einstein-Euler system. We comment further on the 
initial conditions in section 3. 

Our first result is on the vacuum Einstein equations, i.e. Rice = 0: 

Theorem 2.1. Given (E,ho,K) an initial data set, and a relatively compact set U in E, with ho 
in H^f (E) and Kq in i^^(E) with s > | + 2 satisfying the constraints, there exists an Einsteinian 
development M = (U X [0, T],g) satisfying Ricc(g) = weakly. The metric g thus obtained is in 

c°([o,T],H s+1 (u))nc 1 ([o,T],H s (u))nc 2 ([o,T],H s - l (u)). 

Note that the assumption of relative compactness guarantees the uniformity of T. To remove such 
an assumption, either a uniform positivity condition is required on the initial metric or the manifold 
E has to be compact. The proof involves proving existence locally in coordinates and using the finite 
propagation speed to conclude the global result. Details are as in [D]. 

For perfect fluids, the energy-momentum tensor is T = (p + p)u ® u + pg. The Einstein-Euler 
system for a perfect fluid is 

ixlCC^i/ T^^-diiv — J- flu 

(2.1) (A» + p)u a V a u p + {u a u p + g af3 )d a p = 

+ p)V a u a + u a d a p = 
u a u a = — 1 

Notice that we have chosen units such that ^^ = 1. Also note that the first equation maybe written 

as Ricc Mi , = p M „ where p^ u = T^ u - ^p-g^u- 

The Einstein-Euler system also exists for a short period of time: 

Theorem 2.2. Given a smooth function p : K — > (0, oo) with p! > 1, (E, ho,K,po,v) an initial 
data set for the Einstein-Euler system, and a relatively compact set U in E, with ho in H^ (E), 
Kq in iJ^(E), po in iJ^(E), and v in iJ^(E) with s > | + 2 satisfying the constraints, there 
exists an Einsteinian development M = (U x [0, T],g) satisfying the Einstein-Euler system weakly. 
The metric g thus obtained is in C°([0,T], H S+1 (U)) n C l ([0,T], H S (U)) n C 2 ([0,T], H S ~ X (U)), the 
four-velocity u G C°([0,T], H S (U)) n C 1 ([0,T], H s ~ l (U)), and the pressure p e C°([0,T\, H S {U)) n 
C 1 ([0,T],H s ~~ 1 (U)). They obey the initial conditions, hence in particular the orthogonal projection 
of u onto TE is v. 

Remark 2.3. The condition p! > 1 guarantees that the speed of sound is at most that of light. The 
Einstein-Euler system will be a FOSH only as long as this condition holds. This is one advantage of 
having only physical characteristics, as mentioned in the introduction. 

As is well known, these kinds of equations are diffeomorphism invariant and hence an appropriate 
gauge has to be chosen in order to solve them. Traditionally, harmonic coordinates were employed to 
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convert the Einstein equations into a second order hyperbolic system. However, in the Lagrangian 
framework, Choquet-Bruhat and York chose the so-called Cattaneo-Ferrarese (CF) gauge consisting 
of Lagrangian observers following the fluid flow. In other words, a choice of a local orthonormal 
frame e a such that eo = do = u = d x o and the e, = di = A\ (d x j — bjd x o) are Fermi propagated i.e. 
fij = Tqj = 0. Here, local coordinates x a have been chosen. Note that d a should not be confused 
with a coordinate basis (which are denoted by d x a). In this gauge, the Einstein equations were 
re-written to form a reduced FOSH. 

Following [YB] , all the symbols appearing henceforth are to be treated as "abstract" — for exam- 
ple, R a /3[iv is not known, a priori, to be the Riemann tensor of a metric; except, however, for those 
quantities determined at t = 0, in which case they do have their usual meaning. Greek letters refer 
to spacetime coordinates/frames whereas latin ones to spatial coordinates/frames. For the sake of 
brevity, we use the symbol V or a semicolon to mean the "covariant derivative" with the "correct" 
Christoffel symbols i.e. T a Qb = 0, F a bc = u&,r° a = Ya = r a 00 ,T° ab = X ah = Y h aQ . We denote spatial 
covariant derivatives and curvatures with a T and R respectively. A bracket enclosing two letters 
A\ a u indicates antisymmetrisation whereas one enclosing three letters separated by commas (^4[ a ,&,cl) 
indicates cyclic summation A a f, c + At, ca + A ca i,. 

In the case of vacuum, there are no fluid flow lines. Hence we may impose the additional gauge 
choice Y{ = 0. The reduced system of equations as written in [YB] is 











dobi = 


—OjXki 







(2.2) 








Rohij = 9oUJ l h j + Xhk^lj 
RhOiO = —doXhi — XhjXji 
VoRhkX/i = -^kRohX/j, + ^hRokX/i 
Vq-Rq/iA/x = ^iRlhXfi 




where a = A" 1 . 

In the perfect fluid 
reduced system is 


case, let F = J - 


jfpfrp, Poo = l(3p + 


■fi),PiO = 


= 0, and pij - 






doa) = 


—djXhi 














dok = 


-o^Yh 










(2.3) 




d u} % hj + X hk u % kj 
doXfo + XhjXji 


+ Y { X hj 
- Y h Yi - 


— YjXhi = Rohij 

- ViY h - ^-(X hi - 


Xih) = 


-RhOiO 






»'d Y h 


- ViX hl 


-Yi{Xu 


- X a ) + fj,'(Y h doF - 


- x hl d t F) + d hf i'd F 



8ij\(n - p). The 





^oRhkXfi = —VkRohXfj, + V/j-RofcA/x 
^oRohXfi = ^iRihXfj, + Vj/9oi - VoPji 

d fi = -(n + p)X a 

In both vacuum and fluid cases, the constraints are obtained by restricting the Einstein equations 
to t = 0. This enables us to solve for temporal derivatives of all the quantities at t = 0. Substitut- 
ing these expressions in the so-called quasi-constraints 3 which we write below, gives us the actual 



This procedure is necessary because di contains d x o. 
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constraints. The initial connection is the Levi-Civita one. The other quasi-constraints are 

^[h, R i,j]X/i = 
VhRhoxn = Va/v> — V m pao 

Rhkij = Rhkij + XfaXjh — XjfrXih 
— RkhOj = V[ k X h ]j — YjX^h] 

3. Initial data 

The initial data required for the well-posedness of the reduced system is defined on E x {0}: 

• A field of coframes a*- and of covectors bi creating a metric ho on E via h J = aja — Wb 
which is assumed to be positive definite. The initial metric on the manifold M is g(t = 0) = 
-(6 ) 2 + XX^) 2 where 9 i = a)dx^ and 9° = dx° + hdx 1 . 

• Fields ojfj , Xij , and Yj (which is assumed to be zero in the vacuum case) . These are supposed 
to define the Christoffel symbols of the Levi-Civita connection of g initially (with /y = 0). 

• Tensor components Rijtu Roijh an d RoiOj that define the Riemann curvature tensor initially. 

• In the case of the perfect fluid, we also need ji(p) > obeying // > 1, and p > 0. 

In addition, the Einstein equations are imposed on this initial data at t = in order to derive the 
relation between the usual Eulerian initial data of the 3 + 1 decomposition -- which is given in 
theorems 2.2 and 2.1 — and the initial data needed for the FOSH systems, as we now explain. A 
detailed account of the correspondence between initial data sets for the Einstein-Euler system and 
those of reduced equation in Lagrangian coordinates can be obtained by an argument similar to that 
of[D]. 

Given a Riemannian 3-fold (£, ho), choose local coordinates x l on it. Embed it into M = S x 1 as 
E x {0}. Let ej be an orthonormal frame on £ and let 6 l be the dual coframe. Then ho = ^ 9 l ® 9 l . 
In the vacuum case, we may simply define the initial Lorentz metric on S x {0} as g = — (dx ) 2 + h$. 
This corresponds to bi being zero initially. 

In the case of a perfect fluid, define a metric g = hQ + Vidx°9 l — (dx ) 2 on TM restricted to E x {0} 
with Vi being the components of the dual (with respect to ho) of v . This is a Lorentzian metric with 
eo = d x o being a unit timelike vector projecting to v and restricting to ho on E x {0}. Complete eo 
to an orthonormal basis e a . This gives us a*- and bi lying in H S (T,). Calculations similar to the ones 
in [D] maybe used to define co,X, etc on E x {0}. 

The above reasoning combined with the fact that both (reduced) systems (2.3) and (2.2) are 
quasilinear FOSH having initial data in (at least) H s ~ l implies that both systems have solutions in 
C°([0, Ti],H s ~ l (U)) n C^flO, Ti], H s ~ 2 ([/)). The mismatch between the regularity of the initial data 
and that of the solution is then corrected by a bootstrap argument as in [D] following the results of 
[FM]. 

Thus we have a solution to both systems with a) G C°([0,2i], H S+1 (U)) n C^fltyTi], H S (U)) n 
C a ([0,Ti],ir- 1 (Z7)) J b h 4, Xij, Y u p in C°([0,T 1 ],H'(U))nC 1 ([Q,T 1 ],H a - 1 (U)), and R a ^ u in 

c°([o,r 1 ],^- 1 (c/))nc 1 ([o,T 1 ],^- 2 ([/)). 

4. Proofs 

We prove that the constraints and the gauge are preserved. This implies that if the Einstein 
equations are satisfied initially, then they are satisfied in the future. We accomplish these steps by 
proving that the relevant quantities satisfy FOSH systems with zero as their unique solution. Note 
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that by definition, do = -^ and di = A\(-£j — h-^)- We also note that if a linear symmetry of the 
RaPfiv is satisfied initially, then its spatial derivatives are zero. Since the temporal derivatives are 
related to the spatial ones by the evolution equations (which are imposed on the variables at t = 0), 
we see that di applied to such a symmetry also yields zero. 

4.1. Vacuum. The regularity of the quantities that follow shall be addressed later. For now, we 
assume that everything is smooth. Firstly, we see as to why the vacuum Einstein equations are 
implied by the preservation of the constraints and the gauge: 

Lemma 4.1. If the gauge and the constraints are preserved, then Ricc a /3 is zero in the future, if so 
initially. 

Proof. Notice that fiscal = 2divRicc. The given system implies that (assuming the constraints and 
the gauge) RicCo,^ ;7 = RicCo, 7; ^. Contracting a and /3, we see that cfecal = i.e. seal = because it 
is so, initially. Hence 

Ricc j;o = Ricc o;i 
(4-1) * 

RlCC 0;0 = 2^ Rlcc 0J;« 

i=l 

The leading matrix (Mo) for the equations above is 

I i - Bl -B 2 -B 3 \ 

-Bi 1 

-B 2 1 

\ -B 3 1 ) 

where Bi = —A l -bj. It is positive definite (see lemma 11 in [YB]). 

Riccjo-n = RicCjo-i 
(4.2) 

Ricc i0; o = Riccij-j 

Equations (4.1) and (4.2) form a FOSH system (the leading matrix of equation (4.2) also has positive 
eigenvalues by a similar argument as for equation (4.1)). Hence Ricc Qi g = 0. Note that we treated 
Ricc Qj g and Ricc^Q, as distinct variables. □ 

Now, we prove that the constraints and the gauge are preserved i.e., among other things d a 
forms an orthonormal frame for a metric whose Levi-Civita connection's components are T a ab = X a b, 
T^ a = 0, and T^- = uf-, and whose Riemann curvature tensor is R a p^ v . For further use, we prove 
some symmetries of R a /3^ u '■ 

Lemma 4.2. The following relations are satisfied for some time: 

^afifjbv ^a/Sufi, 

Proof. 

9o(RhkOj + RkhOj) = — Vfci?0M)j + VhRokOj + Vfci?0M)j — Vfr-RofcOj 

= 
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Hence RhkOj = —RkhOj (since it is so, initially). Similarly, Rijki = —Rjtkl- We also have 

do(Rijkl + Rijlk) = —Vj(Roikl + Rank) + ^i{Rojkl + Rojlk) 
do(Roikl + Roilk) = ^p(Rpikl + Rpilk) 

The system above is FOSH (having zero as its unique solution). Indeed, the leading matrix is 



/ 1 





B 2 
-B 1 

V o 





1 


B 3 



-Bx 






1 


B 3 



-E>2 —B\ 
B 3 

B 3 

1 
1 



-B 2 







o \ 
-B 1 
-B 2 




1 / 



Its eigenvalues are 1,1± \J B\ + B\ + B\ with multiplicity 2. They are positive for some time (by 
the assumptions on a and b). This means that Rijki = —Rijik and Roiki = —Rouk- Using these 



symmetries of R a /3fj,u, we see that ^(wf ■ + ur ik ) 



x \Hj+^ii)- Hence w^- 



~ u ik- 



U 



By explicit calculation, we see that d$di — dido = —X^d a and that 

dud,] (2fij + X^do + ^dp 



[h"3\ 



where / and c are two new symbols that satisfy fy = and c^- = wf- ., when t = 0. The evolution 
equations for fij and vf '■ = (?- — u^ ., are obtained by differentiating the above equation and using 
the main evolution system. 

<9o(^)<9p =d %dj] - tfjdodp ~ 23 (/ ij )9o - d {X [itj] )d Q - (2/^- + X [itj] )dl 
+ Ro[i,j] P d P + w^Xqdp 

=%a a i] - xfrdadfl - fydpdo - x a p d a ) - 2d G {f ij )d Q - (2/ y + x^jsg 

- 9o(^[jj])i9o + Ro[i,j} p dp + wj^X^dp 

= % d j} d - d {i XX^d b ) - X^(d j} d a + c? aj d p + (2f aj + X M )do) 

- ^(Opdo - X;d a ) - 2do(fi3)do - (2/y + X M )^ - ^(X^])^ 
+ ^[J,j]pdp+U;fyX l ]<9p 

= (-%(Xj) - X^ ., + 41: + flbfcfl, + ofy-Xj)^ 

+ (- X [i,@faj + "^[aj]) ~~ 2 do{fij) ~ Ro[i,j]0 + X^X.j])^ 
Comparing coefficients we see that 

*(*&) = -%,(^) " ^<i] + <$*Z + *ofci]p + ^4 

3o(2/y) = -X" (2/aj + X^.j]) - A)[i,j]0 + XXj] 



(4.3) 



The system (4.3) is a FOSH system for v and / and hence has a unique solution. If the first Bianchi 
identity and the defining equation of the Riemann tensor are satisfied, then / = v = is a solution, 
and hence as promised, X and u> form the Levi-Civita connection of the metric defined by d a . 
Next, we write evolution equations for the first Bianchi identity. In what follows R\ a pu\v = 
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9oR[h,c,d\a — ~~ ^/i-R[c,0,d]a ~~ ^<2-R[/i,0,c]a: — V c -R[ ,/i,d]a: + ^[/i,-Rd,c]Oa 
^0-R[0,ft,,c]a = ViiJpjh^Ja — VlRhcla + V 'oRhcoa 

The above system is FOSH for the variables R\ a f},u]v Indeed, it is symmetric and the eigenvalues 
of the leading matrix are 1 and 1 ± \/Bf + B\ + B^ (which are positive). We will write evolution 
equations for the other terms in the system (4.4). This will prove that the unique solution to the 
above system is (since it is initially). 

Now, we write the evolution equations for VbRhcba — ^oRhcOa and prove that zero is their only 
solution. To accomplish this, we ought to prove that the lower order terms in these equations 
vanish assuming that all the identities (including the Bianchi identities, the Einstein equations, 
^bRhcba ~ Vo-RhcOa = 0, fij = 0, c^ = ujfj - u^ v etc) hold to order zero. 

Vo(X7bRhcbo) = 

Vo(<9b-R/j cW ) — ^bhRacbO — FfcRhabO ~ F bb Rh C aO ~ ^boRhcba) 

= [<9o, db]Rhcbo — dbV c R hbo + db^hRocW — Xbh^iRicbo 

+ ^bh^cRoabO — OJbh^ aRocbO ~ Xbc^lRhlbO + ^bc^ aRohbO 
~ WfcVhRoabO + W M)V cRohaO ~ ^S,V ' hRocaO + -Xft a V c Rohba 

— Xba^hRocba + {RbOhO + XbkXkh)RocbO ~ {Robah ~ Xbk^ k h)RacbO 
+ (RbOcO + XbkXkcjRhoboiRobac ~ XbkU kc )RhabO 

— (Robab - Xbk^>kb)RhcaO + (RbOaO + XbkX ka )R hcba 

= — d c VbRohbO + dh^bRocbO ~ X ba V a RhcbO ~ X b h^lRlcbO 
+ ^bh^cRoabO — ^bhS aRocbO ~ Xbc^lRhlbO + <^bc^aRohbO 
~ WfcVhRoabO + W M)V cRohaO ~ ^S,V ' hRocaO + XbaS? cRohba 

— XbaVhRscba — dbV c RohbO + &bV hRocbO + <9 c Vfei?o/i60 ~ dhNbRocbO 

At this point we note that 

d c VbRohbo - dbV c Rohbo = v^ b d a Rohbo + ^fcb^iRihm + (X c b — Xb c )ViRihbo 

+ ( w cft ~~ U% c )d a RohbO — d c (Xb a RahbO + ^bhRoabO + WfoRohaO + XbaRohba) 
+db{X ca R a hbo + ^>chRoabO + ^cbRohaO + -^ca-Ro/iba) 

Noticing that d b T« k - d c T% k = Rg k + T« k {T p bc - T p cb ) - T« p T p ck + T« p T p k up to to the zeroeth order by 
assumption, we have 

d c VbRohbo - dbV c Rohbo = 

Vcbd a RohbO + Zfcb^lRlhbO + (X c b ~ Xbc)^ iRlhbO 
+ i^cb ~ ^bcW a R()hbO ~ X b ad c R a hbO ~ U bh d c R 0a bo 
— LU bb d c R()haO — Xb a d c Rohba + X ca dbR a hbO 
+ ^ch^bRoabO + ^cb^bRohaO + X ca dbR()hba 
+ (RahbO + Rohba){ — Robca ~ Xb p UJ^ a + X cp UJ ba ) 
+ Roabo(Rabch ~ ^bp^ch + ^cp^bh ~ Xba,X c h + X ca X bh ) 
+ Rohao(Rabcb ~ ^bp^cb ~ ^cp^bb ~ Xba,X c b ~ X ca Xbb) 
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Inserting (4.6) and another equation (the same one as (4.6) with h and c interchanged and the 
sign flipped) into equation (4.5) we see that the zeroeth and the first order terms cancel assuming 
all the identities hold to order zero. 

The other evolution equations are similar. We write only the highest order terms here. The lower 
order ones (indicated by Li) vanish if we assume (as before) that the identities hold to order zero. 

Vo(Vbi? /iOb) = _ diVbRihbo + L\ 
Vo(\7bRhcbk — Vo-Rhcofc) =d c (ViRihok + ^iRmik) 

— dh(ViRi c ok - V/iV;i?co/fc) + L 2 
(4.7) V (ViRi hok + ViRhoik) = dhiV^Ro^) - d c {V {i R ohjtk] ) + L 3 

Vo(V[j i -R /ij,fc]) =di(V[ i Ri hjtk ]) + L,i 
Vo(V[ ,-R/ lc fc,i]) = - d c (V[ R ohk j]) + <9/i(V[ ,-Rocfc,z]) + L $ 

^o(^[0, R Ohk,l]) = ~ ^m(^[0,Rhnik,l]) + L 6 

The system (4.7) (along with equation (4.5)) is easily verified to be FOSH with zero as the unique 
solution if zero initially. 

We note that VhRho\fi an d ^[iRj,k]\^ evolve according to 

Vo(V/ji4oA^) = L-j 

(4.8) 

Vo(V[j : -Rj ;fc ] AM ) = L 8 

Finally, we calculate the evolution of B khj = R kh0j + d k X hj - d h X kj - X pj (io p kh - io p hk ) + X kp b? h - - 
X hp oj p k - and Whkj = Rhkij - Rhkij - X ki Xj h + Xj k X hi (i.e. the definitions of the components of the 
Riemann tensor) 

(4.9) doBkh > = ~ x : iBkha + xlBpk3 - xlBph3 + XpiR[oAh]p 

doWhkj — -^ijRyox^i + u ijBhki + X h iW k uj — X k iWhUj 

The system (4.9) is FOSH having zero as its solution. If B and W are zero then R a p^u is the Riemann 
tensor of the metric. 

Proof of theorem 2.1: Even with the low regularity mentioned in the theorem, we are allowed to 
manipulate expressions assuming that are smooth. This is because, under the regularity assumptions, 
if a FOSH is satisfied weakly then it is so strongly (almost everywhere) [FM], [M]. Equations (4.3), 
(4.4), (4.5), (4.7), (4.8), (4.9) form a FOSH. Using the uniqueness theory for the same, we conclude 
that zero is the unique solution (zero in H s ~ 2 is the same as zero throughout because of Sobolev 
embedding) if the variables are zero initially. A calculation shows that they they are zero initially. 
Such a calculation is quite long and will not be presented here, but it is done in essentially the same 
fashion as in [D]. This also holds for the system in lemma 4.1. This proves that we have a solution 
to the vacuum Einstein equations satisfying all the conditions required by theorem 2.1. □ 

4.2. Perfect fluids. Just as before, we write equations for the preservation of the gauge. Indeed, 
we show that the Levi-Civita connection corresponding to the orthonormal frame defined by d a has 
components Tfj = ufj, r 3 Qi = 0, Tq = Yi, and T^ = Xij. 

Calculations similar to the ones in lemma 4.2 show that the same lemma holds for the Einstein- 
Euler system as well. We assume this implicitly in what follows. We define S a p = Ricc a ^ — p a p so 
that the Einstein equations are S a p = 0. 
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Explicit computation shows that [3o,9j] = Xfd a + Yido, and [di,dj] = (2/^ + X^jjjdo + (vf, + 
u? -,)d p where / = = v initially. If we prove that / = = v is preserved, then indeed the 
components of the Levi-Civita connection are as described above. Computations similar to the ones 
leading to the system (4.3) prove that 



(4.10) 



2<9o/ij = --Ro[«,i]0 + 2X [i,j] / V ij Y a ~ 1fa[j X i]a ~ ^\j, Y i] 



If indeed v = f = 0, and the first Bianchi identity holds, then both the equations in system (4.10) 
are satisfied provided Fij = VuYj-i + Xuji— ^- = 0. Let Pi = Yi + diF. We record the following 
calculations for further use 



(4.11) 



j-,]. 



d\k F, 



[k^fy 



V[i, Y j]+ x li,j]d° F 
d [i P j] -2f ij d F-v? j d p F-u;l J] P a 

L\ + djFhk 



where Lj denote lower order terms as before. They vanish when all the identities are satisfied. The 
evolution of F is given by 



(4.12) 



fi'doFij = lower order + du(doPjj) 

= lower order + d^diX^ - d^d^Xu 

= lower order + did^Xju 

= L 2 + diBiji 



We now write the evolution equation of R\ a fi,n] v as before: 



(4.13) 



^oR[h,c,d]a — -VhR{c,0,d\a ~ ^dR[h,0,c]a ~ ^cR[0,h,d]a + ^[hRd,c]0a 
^oR[0,h,c]a = ^lR[l,h,c] — ^{l,Rh,c]la + ^ ' cSha ~ ^hS ca 



We wish to make sure that the Euler equation Yi = —diF is satisfied. The evolution of Pi is computed 
to be 



d Pi 



(4.14) 



1 ; 

7 - 

1 

Jl' l 



^j X ij ~ ^iXkk — Y (Xki — Xik, 



XaPi 



Sio + (— Riccjo + VjXij — VjXfcfc — Y (X^ — Xik)) 



XuPi 



Now, we calculate the evolution of the (quasi-)constraints (remembering that Rhoio = R-hOiO ~ ^~hi 
where Rhoio is the "true" Riemann tensor). Let B kh j = RkhOj + dkX hj - dhX k j - X p j(uj p kh - u? hk ) + 
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XkpWfrj — XhpU^j — Yj(Xkh — Xhk) and Whkij = Rhkij ~ Rhkij — XkiXjh + Xj k Xhi- 

doBkhj = — XijB kh i + X kl B hl j + X h iB ik j + X p jR[ 0kh ] p 
- Y p Whkpj - YjRo[k,h]o + ZfkhdoYj 
+ v lh d p Y i + Y [kfh]j ~ YjF[kM + ^[K F h]j 
(4.15) =L 3 + d j J : h k 

doWhkij =£4 

do^[iRj,k]av = L 5 

do{V h R hoXll - V A /V) + V^pao) = L 6 
Finally, we compute the evolution of S a p 

Vo^b = V a Sob + VhR[0,a,h]b ~ (^ hRfiQab ~ ^aPbO + ^bPao) 

V S 0a = V k Ska + 2(/i + p)P b 

Vo^ao = VaS'oo + V/ l (i?ohOb _ RobOh) 

Vo^oo = Vfc-Sfco 

Notice that (4.10), (4.12), (4.13), (4.14), (4.15), and (4.16) form a FOSH. They have zero as their 
unique solution if the Bianchi identities, the constraints, and the Euler equations hold. 

Remark 4.3. Notice that the preservation of gauge and the satisfaction of Einstein's equations 
maybe proven in the vacuum case in a manner similar to that of the fluid case. Actually, the vacuum 
case maybe deduced from the fluid case by keeping track of Newton's constant G in the equations 
and setting G = 0. However, we chose to do it otherwise in order to separate the preservation of the 
gauge and constraints from the Einstein equations themselves. This maybe useful in other contexts. 

Proof of theorem 2.2: As before, the regularity does not hinder the manipulations made above. By 
arguments similar to those used in the proof of theorem 2.1, we obtain a solution of the Einstein- 
Euler system. This satisfies almost all the conditions required by theorem 2.2 except ostensibly, 
the regularity of g, because of the regularity of bi that is lower than desired. This however, is an 
artifact of our chosen coordinate system (which depends on v which in turn has lower regularity than 
ho). To see this, first we note that the g we obtained is in C°([0,T],H s (U)) D (7 1 ([0,r],H 5 - 1 (t/')) n 
C 2 ([0, T], H S ~ 2 (U)). However, its restriction to S x {0} is in H s+1 . Using the ADM decompositon of 
lapse and shift, one may choose the initial lapse to be 1, the initial shift vector to be and their time 
derivative appropriately so as to satisfy the wave gauge condition initially. The energy-momentum 
tensor will be in C ([0, T], H S (U)) n C 1 ([0, T], -£P -1 (C/)). The Einstein equations form a quasilinear 
hyperbolic system in the wave gauge. Hence we get a metric solving the full Einstein equations with 
the correct regularity. This coincides our original solution due to local geometric uniqueness [B] . □ 
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